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Charmonium Spectroscopy From Lattice NRQCD 

NRQCD Collaboration presented by A.J. Lidsey ^ 

^'Department of Physics and Astronomy, University of Glasgow, 
Glasgow, G12 8QQ, Scotland, UK 

We present the first set of results for Charmonium Spectroscopy using Non-Relativistic QCD (NRQCD). For 
the NRQCD action the leading order spin-dependent and next to leading order spin-independent interactions have 
been included. Clear signals for the s and p hyperfine splittings have been observed as well as various orbital 
states. 



1. Introduction 

NRQCD is an effective field theory where a cut 
off of the order of the quark mass is introduced 
into QCD in order to remove the mass scale from 
the calculation. Both numerical and theoretical 
work has been done in developing NRQCD on the 
lattice where the lattice spacing provides the cut 
off. The NRQCD action is essentially an expan- 
sion of the usual QCD action in powers of the 

2 

individual quark velocities squared % , where for 

2 

the bottom and charm quarks ~ 0.1 and 0.3 
respectively. For our action we use the lagrangian 
defined in Euclidean space by 

C-NRQCD = -^^Dt^f + ijft__^ 



2Mq 

SCsi + SCsD 

where 

Csi - ~C2i^-^\D.E - E.D)^ + 
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are the next to leading order spin-independent 
terms, and 

CsD = C3^^^n.{D^E - E^D) + 
Q 
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are the leading order spin-dependent terms. Since 
we have included next to leading order spin- 
independent terms it is possible when using 



charm quarks to obtain an accuracy of w 10% for 
quantites such as the spin averaged S-P splitting 
or the 2S-1S splitting. For spin-dependent quan- 
tities like ^Si—^ So an accuracy of only w 30% is 
possible. 

To define Cnrqcd on the lattice Q covariant 
derivatives are represented by shift operators 



and the laplacian is defined by 
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To evaluate the chromomagnetic and electric 
fields the usual cloverleaf expression is used to 
find ¥^^{x) and hence E(x) and B(x). The ar- 
bitary constants appearing in eq. (||) and (||) are 
evaluated by matching NRQCD to full QCD at 
low energies. At tree level all the c's are unity 
1^. The biggest correction coming from loop di- 
agrams will be tadpole corrections arising due to 
the non-linear relation between the gluon field 
and the lattice gluon operator. In order to re- 
move this we use the method suggested by 
where all the U's are redefined by 
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with Uo 0\^TrUpLAQ\0 Since the clover- 
leaf expression involves the evaluation of a pla- 
quette this renormalization will have the effect of 
redefining E and B via 
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which will strongly affect spin-dependent quan- 
tities. Not removing these tadpole contributions 
will result in an underestimation of the electric 
and magnectic fields and hence hyperfine split- 
tings. The tadpoles are the dominant affect in 
radiative corrections so removing them enables 
us to use the tree level values for the c's. 
As with all lattice simulations lattice spacing er- 
rors will need to be investigated and corrected 
for. In NRQCD for heavy-heavy quark systems 
correction terms can be added as a power series in 
the typical quark velocity in the same way as rel- 
ativistic corrections were in eq. and eq. (|^). 
In order to obtain the desired accuracy quoted 
above it is necessary to remove 0{a^) corrections 
from the laplacian by redefining as 

A2=^A+Ar-^^[A+Arp (8) 

i i 

and for the temporal correction it is sufficient to 
redefine the kinetic operator as 

2M 2M 4n AAP ^ ' 

which still allows the evolution equation to be an 
initial value problem, n is an even number in- 
troduced to prevent instabilities occurring in the 
simulation at low quark mass [|l],^ . 

2. Simulation results 

In the simulation we use a lattice of size 12'^x24 
at 13 = 5.7 with quenched gauge configurations 
supplied by the UKQCD collaboration, fixed to 
Coulomb gauge. The operators are smeared using 
a coulombic wavefunction as a smearing function. 
Due to the relatively small size of charmonium it 
is possible to use more than one starting site on 
the lattice and since our evolution equation is a 
simple difference one, different starting times can 
be used too. For our particular lattice we use 8 
different spatial origins and 2 different starting 
times. 

For NRQCD the dispersion relation has the form 
Ev^Mi + ^ + .... (10) 



Splitting 


Lattice result 
(MeV) 


Exp 
(MeV) 




96 (4) 


117 


'P2 -' Pi 


60 (30) 


46 


'Pi -' Po 


52 (9) 


95 



Table 1 

Results for s and p hyperfine splittings for the Is 
and Ip states. 



where, because we have removed the Mg^'tvE' 
term from the original Lagrangian, Mi ^ M2- 
Ml is a redundent zero of energy and it is M2 
which is the important dynamical mass deter- 
mining spin averaged splittings. The bare quark 
mass Mq appearing in the Lagrangian is tuned 
in the simulation so that M2 appearing in the 
dispersion relation for the ^Sq state is equal to 
the experimental value of the r]c. This is done 
by Fourier transforming the correlation function 
into momentum space from which Ep^.^ can be 
found, where Pmin is the lowest allowed non-zero 
momentum on the lattice. Using eq. ( p^ M2 is 
evaluated from Ep^.^ — Eq. 

In our simulation we have used 100 configura- 
tions and a bare quark mass in lattice units of 
0.8. Extracting the lattice spacing from the spin 
averaged S-P splitting, a value known to be in- 
sensitive to the bare quark mass Mq, a value of 
— 1.23 (7) GeV is obtained. Using this we 
find M2 = 3.00 Gev for the ^Sq to be compared 
to the experimental value of 2.98 GeV. 
Shown in figure ([T]) are values for various spin- 
independent splittings where mostly single expo- 
nential fits have been used. The ^D2 — ^ Sq 
splitting has to be compared to the experimen- 
tal ■^Di So result since the ^D2 has not been 
observed experimentally. The 'Di has the same 
quantum numbers as the '^S'l and will appear as a 
third excited state state in that channel making 
it difficult to extract a value for its mass. The 
simulation result is slightly higher as one would 
expect and certainly in the right area considering 
the ^D2 is near the threshold for decay into D 
mesons. For the 2S-1S splitting the 2S was ob- 
tain by performing a three exponential fit to two 
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Figure 1. Results for spin- averaged splittings 
plotted relative to the spin-averaged IS state. 
The vertical scale is in GeV. The horizontal lines 
mark experimental results, the squares are results 
from this simulation. Crosses represent splittings 
which were fixed to experiment. The ■0(3770) is 
a. ^Di state compared to our result for the ^Z?2 
(not spin-averaged). 



correlation functions. The result matches the ex- 
periment although the error bar is large. 
Tabel (|l|) shows values for spin-dependent split- 
tings. For the S hyperfinc splitting the value 
is well within the expected systematic accuracy 
even though it is smaller than the experimental 
value. Since this splitting is caused by a local- 
local interaction the reduction might be caused 
by quenching. This affect has been estimated to 
cause a reduction of « 20% Q. The fact that 
the s hyperfine is within the expected accuracy 
strongly suggests that the improved tadpole cou- 
pling constant in the a.B interaction is the correct 
value to use. For the p hyperfines the statistical 
errors have increased but still a signal can be seen. 
The agreement with experiment suggests that the 
a.D^E term can correctly account for the p hy- 



perfine splitting, again using a tadpole improved 
coupling. 

3. Conclusion 

Using NRQCD it is possible to obtain a spec- 
trum for Charmonium in good agreement with 
experiment. 

In the future more configurations will be used 
to reduce statistical errors. This will provide a 
stringent test of our perturbatively improved ac- 
tion and so of QCD itself. Using perturbation 
theory with the numerical results will then en- 
able a value for the charm quark mass to be de- 
termined. 
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